Light Bending as a Probe of the Nature of Dark Energy 
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We study the bending of light for static spherically symmetric (SSS) space-times which include 
a dark energy contribution. Geometric dark energy models generically predict a correction to the 
Einstein angle written in terms of the distance to the closest approach, whereas a cosmological 
constant A does not. While dark energy is associated with a repulsive force in cosmological context, 
its effect on null geodesies in SSS space-times can be attractive as for the Newtonian term. This 
dark energy contribution may be not negligible with respect to the Einstein prediction in lensing 
involving clusters of galaxies. Strong lensing may therefore be useful to distinguish A from other 
dark energy models. 

PACS numbers: 95.36.+x,95.30.Sf,04.50.+h 



I. INTRODUCTION 

It is still unclear what drives the universe into accel- 
eration recently. While a cosmological constant A is the 
simplest explanation, its value seems completely at odd 
with the naive estimate of the vacuum energy due to 
quantum effects. 

An alternative to A is obtained considering a dynami- 
cal degree of freedom added to the primordial soup. Since 
dynamical dark energy (dDE) varies in time and space, 
its fluctuations are potentially important in order to dis- 
tinguish it from A Another possibility for the expla- 
nation of the present acceleration of the universe is given 
by the geometry itself, through a modification of Ein- 
stein gravity at large distances @: these are known as 
geometric dark energy (gDE) models. A non vanishing 
mass for the graviton is among these possibilities [3j. 

Cosmological observations, such as those coming 
from Supernovae, cosmic microwave background (CMB) 
anisotropies and large scale structure (LSS), have not 
been able to discriminate among dDE models yet (see 
[3| for updated constraints and forecasts). It is therefore 
important to explore observational tests at astrophysical 
level for objects which are detached from the cosmologi- 
cal expansion. 



II. DEFLECTION OF LIGHT 

We shall consider a SSS metric in terms of the physical 
radius r 

ds 2 = -B{r)dt 2 + A(r)dr 2 + r 2 dti 2 . (1) 
Such a metric with 

B(r) = A- 1 (r) = l -r 2 , 2 

r 6 

describes the Schwarzschild-de Sitter (SdS) space-time, 
the vacuum solution of Einstein equations in presence of 
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FIG. 1: Deflection of Light 

a cosmological constant A (in c = 1 units, where c is 
the velocity of light). With the above metric, the classic 
general relativistic tests can be computed in analogy with 
the Schwarzschild (henceforth S) textbook case |f|. We 
shall restrict here to light bending, leaving other results 
for elsewhere 0, 0] ■ 

The textbook deflection angle for a photon in a SSS 
metric can be easily extended by keeping into account 
the finite distance r at which the object is located as (see 

Figm 



<p = 2|^(r ) - 0(r)| 

= 2 1(r, tq) — 7r + 2 arcsin 
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with ro representing the minimal distance between the 
geodesic of the light and the lens [1, Q and I(r, ro) given 
by: 



r ^V2( r ') 



B(r ) 
B{r') 



- 1 



-1/2 

dr' 



l(r',r )dr'. (4) 
In the textbook treatment Q , both the distances of the 
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observer and the source from the lens are much greater 
than r and the upper limit of integration can be taken as 
oo. For DE metrics with non-asymptotically flat terms, 
it is safe to keep finite the upper limit of integration. 

By inserting the SdS metric in X(r, ro), it is easy 
to check that the terms including A cancel out! This 
cancellation is due to the particular form of the non- 
asymptotically flat terms of the SdS metric and physi- 
cally means that A is truly a potential offset for a massless 
particle (see below). The light bending angle is therefore 
4GM/r even in presence of A [ToL fTlj , although all the 
other tests and kinematical quantities differ from the S 
cased 0- 



non-perturbative solution found by Vainshtein (V) 
for a massive graviton [131 ] with mass m g . Massive 
gravity (MG) is an alternative to a cosmological 
constant [H; 

C = 1 +3GM 7l , D = 2GM + 3G 2 M 2 - n , a = 
1 7 7i =72 =7 corresponds to the general SSS 
in conformal gravity [lj] (see also [15| for linear 
correction to the Newtonian potential) . Conformal 
gravity contains the SdS solution without adding 
any cosmological constant to the Weyl action, but 
includes terms linear in r as well. Considering C ~ 
1 and D ~ 2GM is a good approximation for the 
values of 7 we are interested in; 



III. GEOMETRIC DARK ENERGY STATIC 
SPHERICALLY SYMMETRIC METRICS 

In order to study the contribution of DE to light bend- 
ing, we need physically motivated SSS metrics which dif- 
fer from the SdS case. As an example, we consider SSS 
metrics parametrized by: 



B(r) 
A-\r) 



r 

C - - - 7l r° 
r 



(5) 



This paramctrization covers SSS metrics of well studied 
gDE models: 



C 



1, D 



2GM , 



3/2, 71 



372/2 = -m 2 gV / 2GM/13 [12| corresponds to the 



C = 1, D = 2GM 
y/GMjrJ, corresponds 



* = 1/2,7! = 72/2 = 
to the self-accelerating 
of the brane induced _gravity Dvali- 
1711. r c is the 



branch [l6| 

Gabadadze-Porrati (DGP) model 
crossover scale beyond which gravity becomes five- 
dimensional. 



In the parametrization (T5J) we have omitted a possible 
r 2 term, since we have already shown that it does not 
contribute to light bending. For a = 1 the SSS metric is 
valid till to the particle horizon radius rjj ~ 7 _1 (for the 
SdS metric this radius is (3/A) 1 / 2 ). In the other two cases 
the metric is roughly limited by the so-called V radius 
ry , which denotes the scale at which the Newtonian term 
becomes comparable to the non-asymptotically flat term, 
i.e. r v ~ (GAf/ 7 ) 1/(1+Q) with 7- 1 = min {7^ , 72 -1 }. 
At this scale we expect deviations from Einstein gravity. 



J 



We now expand the integrand considering the DE terms in Eq. ((4]) as second order with respect to the Newtonian 
correction, and, defining /(r, ro) = ro) + Ide (r, ro) we obtain: 

1-1/2 



W,r ) = / — 

J rn r 



-1-1 

ro 



1 



GM 



r' ro(r' + ro) 



GMr' 3(GM) %{GM)" 



2r' 2 



3(GAf) 2 r' 

ro(r' + r ) 2rg(r' + r ) 



^DE(r,r ) = / ^7- 



1-1/2 



-) -1 

ro 



[(71 - 7 2 )r'( Q+2 ) + 7 2 r /2 r£ - 71 r' a r { 
2(r' 2 - r 2 ) 



It is possible to analitically compute both the terms in Eqs. (|6l7p : 
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(6) 
(7) 

(8) 



with 



i a (y) 



/DE(r,r ) = f I a (L) ( 71 - a72)+7 2 f \ 
2 VoJ 2 ^(r/ro) 2 - 1 



dx- 



2 r(i - f 



l 2Fl 2 



1 1 — a 3 — a 1 



(9) 



(10) 



It is important to note that the physical structure of three values a. Of course, the type of DE contribution 
gDE metric kill the first term in Eq. ([9]) for all the 
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to the deflection angle is not of the parametrized post- 
Newtonian (PPN) form [|], since the correction in the 
metric coefficient is not of the PPN form. Note that our 
non- vanishing result for the DGP model corrects previous 
claims in the literature [18l | and agrees with previous re- 
sult for Weyl gravity when r — > oo [l9l | . A non-negligible 
DE contribution to the total deflection angle is expected 
for large virialized objects, i.e in the case of clusters of 
galaxies. 

In the following, we split the Einstein term ipE from 
the dark energy correction Aip to the total bending angle 
ip. Therefore, the DE contribution has to be compared 
with the S term assuming r — ► oo (up to the second order 

My- 

<~Pe = + V {U) = 



This comparison is shown in Table [J 

The sign of the DE contribution can be understood by 
the one dimensional motion for the photon: 



with dX = r dtf/J and 



Kff(r) 



2r 2 



Mr) 



(7i ~ l2)r a ' 
J 2 B(r) 



2r 2 



C 



D 



7l r- 



(71 - l2)r a 



J 2 (C 



D 



l2 r a ) 



This potential leads to the following force F on a photon: 



F(r) 
J 2 



off 



J 2 dr 



C 



3D 
2^ 



+ ^(a~2)r Q - 3 + ( 72 - 7l ) 



2 J 2 (C-2 



72 r" 



C D 



(13) 



From the above it is again clear how A acts as a null force on the photon (in agreement with the null contribution 
in light bending [2l|). When a < 2, in addition to the standard terms, the third term acts as an attractive force 
for 71 > 0. We also note that the sign of the fourth term is model dependent 22]. Although DE in cosmology is 
associated to a repulsive force (which accelerates the expansion) , in a static configuration it may add to the Newtonian 
mass term in deflecting light, explaining the positive contribution which we find in Eq. 











A<^Wcyl 


A(^DGP 


Galaxy (10 11 M ) 


2.0 x icr 6 


1.9 x icr 12 


— 4.6 x icr 11 


1.0 x 10~ 6 


3.9 x 10~ 9 


Galaxy Group (10 13 M Q ) 


1.9 x 10" 5 


1.9 x icr 10 


-7.2 x 1(T 9 


1.0 x 10~ 5 


2.3 x 10~ 7 


Cluster (10 15 M Q ) 


1.8 x 10" 4 


1.8 x 10" 8 


-1.1 x 10" 6 


1.0 x 10~ 4 


6.6 x 10" 6 



TABLE I: Expected deviations of the three theories considered from Einstein theory. The parameters used are the following: 
m g — 10~ 31 eV, 7 _1 = 10 Gpc, r c — 5 Gpc. The distance of closest approach are ro = 10, 10 2 , 10 3 Kpc for galaxy, galaxy 
groups, clusters, respectively. The upper limits of integration have been taken as rv for a = 1/2 ,3/2 and 00 for Weyl gravity. 
We remind that the Hubble distance 1/Hq is ~ 4 Gpc for Ho = 72 km s _1 Mpc -1 . As written also in the text, the contribution 
Ai^bv coming from the region beyond ry is taken into account only in DGP model [l6l . [l8j] and it is given by Ai^bv = 
{G{r max ) - G(r v )) 2GM/r where G(r) = (l - (r- /r) 2 ) 1/2 (1 + /(w)(r /r)) (1 + (m/r))" 1 with f(u) = (w+l)/(2w + 3). The 
parameter co = — 3r c Ho is set to —15/4 and r ma x = 1 Gpc. 



We note that the validity of the SSS metrics for a — 
1/2,3/2 up to the V radius - and not up to the particle 
horizon - can be an important limitation to the appli- 
cability of our findings. In the DGP (massive gravity) 
model with r c — 5 Gpc (m g ~ 10~ 32 eV), the V radius 
for the Sun is 3.2 x 10 18 m (7.5 x 10 20 m), and therefore 
much larger than the size of the solar system ~ 6 x 10 12 



m (taken as the size of the Pluto orbit). For clusters 
with mass ~ 10 15 M Q , instead, the V radius is 10 Mpc 
for DGP and 24 Mpc for MG (with the same parame- 
ters used above): such radii are remarkably close to the 
intercluster distance, i.e. O(10)Mpc. It is then clear 
that an understanding of the SSS metrics beyond the V 
scale (the so-called matching problem [l2|) is needed for a 
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of the source measured respect to the line from the ob- 
server to the lens; d^s and dos are the angular diameter 
distances between the lens and the source and between 
the observer and the source, respectively. On considering 
for simplicity alignment between the lens and the source, 
an Einstein ring forms with angle Be — B{j3 = 0). From 
our results, Be is affected by both the non-perturbative 
SSS potential around the lens (ip ^ AGM /r if DE ± A) 
and the cosmology of a given model. The gDE correc- 
tions to the Einstein deflection angle for clusters in Eq. 
© are as important as the cosmology for an observable 
as Be- The differential of Be is 



FIG. 2: Dark Energy contribution A</?de to the bending angle 
(in radiants) as a function of GM/ro- The solid line (A</jde = 
0) is the Einstein prediction for A; the dotted, long-dashed 
and dashed are the DGP, Weyl and MG models, respectively, 
with ro = 1 Mpc (the other parameters are those of Table [I}. 
Note that 7 is M independent [14| . 

quantitatively exact calculation of light bending in these 
models. While the MG metric is not known beyond ry, 
the DGP metric beyond ry admits a scalar-tensor de- 
scription [l|| [3 , which we use for the values reported 
in Table U (see the caption) and in Fig f5] where we show 
the prediction of the considered models compared with 
Einstein theory as function of GM/ro. 

IV. DISCUSSIONS AND CONCLUSIONS 

We have discussed light bending in SSS with DE. The 
importance of general relativistic tests, such as the per- 
ihelion precession, has been already emphasized for the 
DGP model [16, 23], while little attention was previously 
paid to light bending. These tests are complementary to 
the observational signatures of dark energy in cosmolog- 
ical context, mainly based on the behaviour of perturba- 
tions. In cosmology, an important difference between A 
and dDE or gDE is the presence of DE perturbations in 
the latter case, which are at least gravitationally coupled 
to the other types of matter. Such DE perturbations are 
therefore a key point to distinguish A from dDE or gDE 
in CMB and LSS, and sometimes may become so impor- 
tant to strongly constrain models [j], 0, [25| with respect 
to what Supernovae data can do. 

In this article we have shown that in objects which have 
detached from the expansion of the universe, A may be 
distinguishable from other DE models through the bend- 
ing of light. In order to link our findings with observa- 
tions, we should insert <p in the lens equation, e.g. [2dT |: 

e-0 = ^ (w) 

dos 

where B and (3 are the angular positions of the image and 



^£ = ^ + aln £i, (15 ) 
Be <P dos 

which reveals how cosmological information is encoded 
just in the second term to the right. By considering the 
cosmology of the DGP model for instance [18j |, one finds 
that the second term is ~ — 0.06(A12m/^m) + AH /H 
for a source and a lens located at z = 1 and z = 0.3, 
respectively (fijyi ~ 0.3 and Hq are the present matter 
density and Hubble parameter, respectively). On con- 
sidering the uncertainties on the cosmological parameters 
of the order of percent, this simple quantitative example 
shows how the corrections to the Einstein deflection an- 
gle we have found in Table I should be taken into account 
in the study of strong lensing by clusters. 

We believe that results similar to what we have found 
here for gDE models, might occur for dDE scenarios as 
well, in which the non-asymptotically flat term is due 
to the non-perturbative clumping of DE into objects de- 
tached from the cosmological expansion. However, dDE 
models may be less predictive than gDE models: gDE 
contain the same number of parameters of ACDM, while 
dDE may need more. Let us end on noting that some of 
the gDE models considered here may have serious theo- 
retical issues [I?], HH , whose resolution clearly go beyond 
the present project. However, the main result in Eq. 
© of this paper remains valid: models alternative to 
general relativity with a cosmological constant predict a 
correction to the Einstein angle, which can be used to 
distinguish A from other DE models. 
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